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A self-consistent system of interacting spinor and scalar fields is considered within the 
scope of Bianchi type VI cosmological model filled with a perfect fluid. The contribution 
of the cosmological constant (A-term) is taken into account as well. Exact self-consistent 
solutions to the field equations are obtained for a special choice of spatial inhomogeneity 
and the interaction terms of the spinor and scalar fields. It has been found that some special 
choice of metiic functions can give rise to a singularity-free solutions independent of the 
value and sign of the A term. It is also shown that the introduction of a positive A, the 
most widespread kind of dark energy, leads to the rapid growth of the universe, while the 
negative one, corresponding to an additional gravitational energy gives rise to an oscillatory 
or non-periodic mode of expansion. The role of the spatial inhomogeneity in the evolution 
of the universe is clarified within the scope of the considered models. 
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I. INTRODUCTION 



The role of the spinor field in the evolution of the Universe has been studied by a number of 
authors 01 0, ll, III 1, [(I 0]. A principal goal of those studies was to find out singularity-free 
solutions to the corresponding field equations. The gravitational field in these cases was given by 
an anisotropic Bianchi type-I (BI) cosmological model. It has been found that the introduction of 
the nonlinear spinor field in some special cases results in a rapid growth of the Universe. In light of 
that study many authors believe it is possible to consider the spinor field as one of the candidates 
to explain the late time accelerated mode of expansion ll§.[^.floll. 

Though as an anisotropic cosmological model cosmologists basically consider Bianchi type-I 
space-time, there are still a few other models that describe an anisotropic space-time and generate 
particular interest among physicists lllll. Il2l OS Q H M 0, Hi. In Ref. d methods of 
dynamical systems analysis were used to show that the presence of a magnetic field orthogonal 
to the two commuting Killing vector fields in any spatially homogeneous Bianchi type VIq vac- 
uum solution to Einstein's equation changes the evolution towards the singularity from collapse to 
bounce. The authors in Ref. [13] studied the problem of isotropization of scalar field Bianchi mod- 
els with an exponential potential(s). Other papers mentioned above are devoted to tilted perfect 
fluid solutions, chaotic singularities, and conditional symmetries. 
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In a recent paper 111911 we studied the self-consistent system of the nonlinear spinor field and 
an anisotropic inhomogeneous gravitational field in order to clarify the role of the spinor field 
nonlinearity and the space-time inhomogeneity in the formation of a singularity-free universe. As 
an anisotropic space-time we chose a Bianchi type- VI (BVI) model, since a suitable choice of its 
parameters yields a few other Bianchi models including BI and FRW universes. It can be noted 
that unlike the BI universe, the BVI space-time is inhomogeneous. Inclusion of inhomogeneity 
in the gravitational field significantly complicates the search for an exact solution to the system. 
The purpose of this paper is to study a self consistent system of spinor, scalar and Bianchi type- VI 
gravitational fields in presence of a perfect fluid and a cosmological constant and study the role of 
corresponding material fields in the evolution of the universe. 

II. BASIC EQUATIONS AND THEIR GENERAL SOLUTIONS 

We shall investigate a self-consistent system of nonlinear spinor and Einstein gravitational 
fields. These two fields are to be determined by the following action: 



^(g;i/^) = J &yf=gd£l (2.1) 

with 

— =2g ~t~ =^sp ~{~ "^sc ~&'vaX "t - =^pf- (2.2) 

The gravitational part of the Lagrangian (12.21 ) Jzf g is given by a Bianchi type- VI (BVI hereafter) 
space-time, whereas the terms Jzf sp , Jz? sc , and Jzfj nt describe the spinor and scalar field Lagrangian 
and an interaction between them, respectively. The term Jzf p f describes the Lagrangian density of 
the perfect fluid which minimally couples to the spinor and scalar fields through gravitational one. 

A. Matter field Lagrangian 

For a spinor field i//, the symmetry between y and \fr appears to demand that one should choose 
the symmetrized Lagrangian [20]. Keeping this in mind we choose the spinor field Lagrangian as 



JSf, -- 



Myfryr. (2.3) 



Here M is the spinor mass, is the covariant derivatives acting on a spinor field as 112 ll 12211 



where are the Fock-Ivanenko spinor connection coefficients defined by 

r n = \f( r hrv-dn7o)- (2-5) 
The massless scalar field Lagrangian is chosen to be 

&*c = \<P,aV a - (2.6) 
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The interaction between the spinor and scalar fields is given by the Lagrangian ^\ 

^m = \<P,aV a F. (2.7) 

Here A is the self-coupling constant and F = F(I,J) is some arbitrary functions of invariants I = S 2 
and J = P 2 generated from the real bilinear forms S = ty/ty/ and P = ixfry 5 i/a of the spinor field. 

The contribution of the perfect fluid to the system is performed by means of its energy- 
momentum tensor, which acts as one of the sources of the corresponding gravitational field equa- 
tions. So here we do not need to write the Lagrangian density Jz? p f explicitly. The reason for 
writing Jf p f in Eqs. (12.11 ) and (12.21) is to underline that we are dealing with a self-consistent sys- 
tem. An interesting discussion on the action and Lagrangian for a perfect fluid can be found in 
Refs. EHm. 



B. The gravitational field 



=£^grav — ? (2.8) 



The gravitational part of the Lagrangian in (12.21 ) has the form: 

R 

IK 1 

Here R is the scalar curvature and K is Einstein's gravitational constant. The gravitational field in 
our case is given by a BVI metric: 

d s 2 = dt 2 - a 2 e- 2mz dx 2 - b 2 e 2nz dy 2 - c 2 dz 2 , (2.9) 

with a, b, c being functions of time only. Here m, n are some arbitrary constants and the velocity 
of light is taken to be unity. It should be emphasized that the BVI metric models a universe that is 
anisotropic and inhomogeneous. A suitable choice of m, n as well as the metric functions a, b, c 
in the BVI metric given by (12.91) generates the following Bianchi-type universes: (i) for m = n 
the BVI metric transforms into a Bianchi-type V (BV) universe; (ii) for n = the BVI metric 
transforms into a Bianchi-type III (Bill) universe; (iii) for m = n = the BVI metric transforms 
into a Bianchi-type I (BI) universe and finally, (iv) for m = n = and an equal scale factor in all 
three directions the BVI metric transforms into a FRW universe. 

The metric (12.91) has the following nontrivial components of Riemann and Ricci tensors: 



p 01 _ d r,02 _ b p03 _ c 
K 01 — K 02~ ~Ti A 03 — - ~> 

a b c 



p l2 _ mn ab 13 _ m 2 cd 23 _ n 2 be 

K 12 — ? 7' ^ 13 — o 5 A 23 — ~ T~ T~- 

c L ab c l ca c 1 be 



r,n ( ® b , . c 

/?3 = m n (m — n)- 

V a b c 



— —[ ^ + To + ^ 



a b c 



R\ = ( t - — 

\a ab ac c z 



a db dc m 2 — mn 



b db be n 2 



R 2 = + + 



mn 



^b ab be c 2 

;c dc be m 2 + n 2 
Ri = - - + — + 



c ac be c 



2 
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We write the components of the Riemann and Ricci tensors as invariant characteristics of space- 
time which one needs to know in order to investigate the existence of a singularity (singular point) 
are composed of these tensors together with the metric tensor. Although in 4D Riemann space 
there are 14 independent invariants HI 26], it is sufficient to study only three of them, namely 
the scalar curvature I\ =R,h = fyv^ MV and the Kretschmann scalar I 3 = R a ^v R<xptlv 10, 0]. 
From the Riemann and Ricci tensors written above one finds 



h 
h 



R = — 



(R° ) 2 +(R\y+(R 2 2 y+(Riy+R° 3 Rl 

4[ (R%) 2 + «) 2 + « 3 3 ) 2 + (R^) 2 + (4\) 2 + {R 2 hY 



x — dbc — abc 
i\ 2 , f D 2\2 



dbc 



ab 



(m 2 — mn + n z 



3\ 2 



(2.10a) 

(2.10b) 
(2.10c) 



where we define 

x = abc. (2.11) 

From (12.101) it follows that I\ oc 1 fx, I 2 l/x 2 , and It,<^ l/x 2 . Note that the remaining 1 1 invariants 
are composed of two or more Ricci and/or Riemann tensors and hence are inversely proportional 
to (x) r , where r is the number of tensors in the corresponding invariant. Thus we see that at 
any space-time point where X = 0, the invariants I1J2J3 become infinity; hence the space-time 
becomes singular at this point. 



C. Field equations 

Let us now write the field equations corresponding to the action (12.11) . 

Variation of Eq. (12.11) with respect to the spinor field y(\fr) gives the following spinor field 
equations: 

*'7 M V ^xj/ - My + Sty + #iyV = 0, (2.12a) 

N^xj/Y^+Mifr-^ifr-^iifry 5 = 0, (2.12b) 

where we use the notation 

^ n dF X fydF n dF X „dF 

9 = XS(p ya (f)> a -^j = — , ^ = XP M' a -jj = 2^ a< P dp' 

Since the nonlinearity in the foregoing equations is generated by the interacting scalar field, Eqs. 
(12.121) can be viewed as spinor field equations with induced nonlinearity. 



Variation of Eq. (12.11) with respect to the scalar field yields the following scalar field equation: 

1 d 



(v / =gg VM (l+A J F)(p, M )=0. (2.13) 



Finally, varying Eq. (12.11) with respect to metric tensor g^ v one finds the Einstein's field equa- 
tions. On account of the A term they have the form 

Rj l -ls;R=KT^ + d;A. (2.14) 
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where Ry is the Ricci tensor, R is the Ricci scalar, and Ty is the energy-momentum tensor of the 
matter fields. In our case, where space-time is given by a BVI metric (12.91 ), the equations for the 
metric functions a, b, c read 



v-t! -1- A 


\Z.. 1 Jd) 


kt£+a, 


(2.15b) 




(2.15c) 


kT§ + A, 


(2.15d) 


KTl 


(2.15e) 



b c be n 2 

U + 2 

b c be c L 
c d cd m 2 

- + - + J" 

c a ca c z 

d b db mn 

a b a b c L 

db be cd m 2 — mn + n 2 

+ + 2 

a b b c c a c z 

d b . . c 

m n (m — n)~ 

a b c 

Here over dots denote differentiation with respect to time (t). The energy-momentum tensor of the 
material field is given by 

T i = T fi (s P ) + T ii (sc) + r M fint) + r /x ( P f) ■ ( 2 • 1 6) 

Here ^J S p) is the energy momentum tensor of the spinor field defined by 

T nU = \s pv [w7^vW+ ¥7v^n¥- VnVTvY- VvWuV) (2-17) 
The term Jz? sp in view of Eq. (12.121) takes the form 

j£? sp = -(0S + SfP). (2.18) 
The energy momentum tensor of the scalar field is given by 

y sc) = W V -^sc. (2.19) 

For the interaction field we find 

^(int) = ^ < P,M < P' V - 5 ;^inf (2-20) 

^M(pf) * s tne ener §y momentum tensor of a perfect fluid. For a Universe filled with a perfect fluid, 
in a comoving system of reference such that u^ 1 = ( 1 , 0, 0, 0) we have 

r M(pf) = (P + e)«/-¥P= (e> "P. "P. -P)- (2-2D 

The energy e and the pressure p of the perfect fluid obey the following equation of state: 

P = £e, (2.22) 

where £ is a constant and lies in the interval £ G [0, 1]. Depending on its numerical value, £ 
describes the following types of Universes ll29ll : (i) £ = (dust Universe); (ii) £ = 1/3 (radiation 
Universe); (iii) £ e (1/3, 1) (hard Universes) and (iv) £ = 1 (Zel'dovich Universe or stiff matter). 



6 



Bijan Saha bvissl.tex February 7, 2008 



Here we note once again that the perfect fluid is minimally coupled to the system. Being one of 
its sources, the perfect fluid leaves its trace on the gravitational field which, in turn, influences the 
behavior of the spinor and scalar fields. 

From (12.51 ) we find the following spinor connections for the metric (12.91) 



r = 0, T, = \t[df-m-fY m \ r 2 = ~f[bf + n~f]e n *, T 3 =\cff. 



It is easy to show that 



1 t _ n m — n a 

2 x 2c 



The Dirac matrices y" (x) of the curved space-time are connected with those of Minkowski 
space-time as follows: 



y° = y°, y l =fe mz /a, f = f/be nz 1 y 3 = f/c, 



with 



7° 



/ 

-/ 



f 



where <7; are the Pauli matrices: 

.i _ / 1 
1 



& 
-& 



-i 

1 



r 5 = r 5 



-/ 
-/ 



1 
-1 



Let us consider the spinors and the scalars to be functions of t only, i.e., 

Under this assumption for the the scalar field from (12.131 ) we find 

clt 



+ XF) 



Co 



const. 



For the spinor field from (12.121 ) we obtain 



r°0+^)-( 



m — n 

2c 
m — n 

2c 



y 3 i// + i^i/z + ^y 5 i/A = 0, 
ty/y 3 — i&y — xfry 5 = 0. 



Here we define 4> = M — 3) . Let us introduce a new function 

Uj(t) = VT\]/j(t). 

Then for the components of the nonlinear spinor field from (12.251) . one obtains 

~m — n 



u\ + iQu\ — 
ii2 + i^U2 + 
U3 — z'4>i/3 — 
t<4 — /4>«4 + 



2c 
m — n 

2c 
m — n 

2c 
m — n 

2c 



u 3 
114 

U\ 

u 2 



0, 
0, 
0, 
0. 



(2.23) 
(2.24) 

(2.25a) 
(2.25b) 



(2.26a) 
(2.26b) 
(2.26c) 
(2.26d) 
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Using the spinor field equations (12.121) and (12.251) . it can be shown that the bilinear spinor forms 

S= yy = vv 1 P = ixfry 5 ^ = iv*fv, A = yry 5 y° i// = vy 5 y°v, 
A 3 = ^ry 5 y 3 i// = i7y 5 y 3 v, V° = xfrfx^r = vy°v, V 3 = xfrf 3 y = v^v, 
Q 30 = ixj/'f'fxj/ = iv'ffv, Q 21 = \frff i y 5 \ir= i\frf 1 f 1 \i/= ivf-fv, 

obey the following system of equations: 



S - 2&Aq 
Po-24>A° 


= o, 
= o, 


(2.27a) 
(2.27b) 


A o ~ ^~^A\ + 2<PP + 20 S 


= o, 


(2.27c) 


. , m — n n 
A o c A 


= o, 


(2.27d) 


yO m-n 3 


= o, 


(2.27e) 


V Q 3 m ~ n V + 2^Ql° 20 (ff 
c 


= o, 


(2.27f) 


Go° -2 W 3 


= o, 
= o, 


(2.27g) 
(2.27h) 



where we use the notation Fq = tF. Combining these equations and taking the first integral one 
gets 

(5 ) 2 + ( J P ) 2 + (A[j) 2 -(A 3 ) 2 = d = const, (2.28a) 
(V 3 ) 2 + (^°) 2 + (^ 1 ) 2 -(y ) 2 = C 2 = const. (2.28b) 

Now let us solve the spinor field equations (12.261) . From the first and the third equations of the 
system (|2.26l) one finds 

"13 = (& - <2)4s - 2^"i3 + (& + £>), (2.29) 

where, we denote uu = u\/u?, and Q — [m — n]/2c. Equation (12.291) is of the Riccati type [30] with 
variable coefficients. Transformation [31] 

vi3 = exp(-yV-G)"i3*), (2.30) 
leads from the general Riccati equation (12.291) to a second order linear one, namely, 

- <2)vi3 + [2/<&(Sf - Q) - & + Q] via + (Sf - 2) 2 (^ + fi)vi 3 = 0. (2.31) 

Sometimes it is easier to solve a linear second order differential equation than a first order nonlinear 
equation. Here we give a general solution to (12.291) . For this purpose we rewrite (|2.29l) in the form 

>vi3 = - G)w? 3 <T 2, ' J '*(' )<ft + (Sf + 2) e 2 ' (2.32) 

where we set M13 = W13 exp[— 2i J <&(t)dt]. This is an inhomogeneous nonlinear differential equa- 
tion for W13. The solution for the homogeneous part of (12.321) . i.e., 

vvi3 = {&- 2)w 2 3 exp f -2i [ <f>(t)dt) (2.33) 
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W13 



J{& - Q) exp f-2i j®{t)dt\ dt + C 



(2.34) 



where C is an arbitrary constant. Then the general solution to the inhomogeneous Eqn. (12.321) can 
be presented as 



W\3 



J(&-Q) exp \-2i J <Z>(t)d?J dt + C(t 



(2.35) 



with the time dependent parameter C(t) to be determined from 



C 



J(&-Q) exp^-2i J <t>(t)dtjdt + C(t) 



' + Q)e 



2if&(t)dt 



(2.36) 



Thus given a concrete nonlinear term in the Lagrangian and the solutions of the Einstein equations, 
one finds the relation between u\ and u?, («2 and U4 as well), hence the components of the spinor 
field. 

Now we study the Einstein equations (12.151) . In doing so, we write the components of the 
energy-momentum tensor, which in our case read 



7/ ° = MS+^(p 2 (l+XF)+£ 1 

t} = r 2 2 = r 3 3 = @s+&p--(p 2 (\+?iF)-p. 

Let us demand the energy-momentum tensor to be conserved, i.e., 



i V.UT 1 «,i J V 1 VU 1 R — u - 



T \i. 



(2.37a) 
(2.37b) 

(2.38) 



Taking into account that is a function of t only and = T 2 = T 3 3 , from (12.381) we find 



T ° + Z [lg-T})=0. 



(2.39) 



In view of the scalar field equation and 4>5o — ^Pq = which follows from (|2.27l) the Eq. (12.391) 
yields 

(2.40) 



(2.41) 



£ = --[£ + p). 



Further using the equation of state (EOS) (12.221) for e and p one finds 



£0 



Ceo 



Let us return to Eqs. (12.151) . In view of (12.371) . from (12.15el) one obtains the following relation 
between the metric functions a, b, c: 



-J =^\-) ) -4^ = const. 



(2.42) 
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Subtracting (I2.15al) from (|2.15b|) we find 

d r d 
dt 



MS)}] 



2 2 

m — n 



-X. 



dt k \ £> / J J c- 
Analogously, subtraction of (I2.15al) from (I2.15cl) and (I2.15bl) from (|2.15cl) gives 



(2.43) 



and 



d r r, /a 
— T— < In - 
dt i dt I \c 



d r d < /b 
dt\ \c 



dt 



mn + n z 



mn + m z 



(2.44) 



(2.45) 



respectively. It can be shown that, in view of (12.271) and (12.421) . the Eqs. (12.431) . (12.441) . and (12.451) 
are interchangeable. 

Taking into account that x = abc, from (12.421) we can write a and b in terms of c, such that 



Jfx n c m - 2n 



and 



-in n—2m 



1 / (m+n) 



l/(m+n) 



(2.46) 



(2.47) 



Thus we find a and £ in terms of c and x. In doing so we employed only four out of five Einstein 
equations, leaving (|2.15d|) unused. Addition of (|2.15a|) . (|2.15b|) . (|2.15c| ). and (|2.15d| ). multiplied 
by 3 gives 



T m 2 — mn-\-n 2 3k, n m1l 

- = 2 ? + T [r ° + 2?] + 3A, 



which in view of (I2.37I) takes the form 

T m 2 — mn + n 2 3fC 



+ ^[M5 + ^5 + ^P+ (1 - C)£o/? 1+? ] +3A. 



(2.48) 



(2.49) 



As one sees, we only have one equation with two unknowns. In order to resolve this problem, 
we have to assume c as a function of x (or vice versa). Given a concrete form of the spinor field 
nonlinearity one finds the solution of (12.491) . This is exactly what we do in the next section. 



IH. ANALYSIS OF THE RESULT 

In the preceding section we derived equations for the spinor, scalar and gravitational fields 
and their general solutions. Comparing the equation with those in a BI universe (see e.g., Ref. 
[Hi@]) we conclude that introduction of inhomogeneity in gravitational (through m and n) imposes 
additional restriction on the metric functions. In fact, Eq. (I2.15el) which connects the metric 
functions a, b, c among themselves, does not figure in the BI universe. In the foregoing sections 
we obtained the solutions to the field equations in terms of x, whereas, the equations for x contains 
the function c explicitly, i.e., we have just one equation for two unknowns x and c. In order to 
resolve this we have to impose some additional condition relating x and c. Though this assumption 
imposes some restrictions on the metric functions, though leaving the space-time anisotropic. 
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Let us consider the case when F = F(I) = F(S) only. In this case from (12.271) we obtain 

S = C /t, (3.1) 

with Co being some arbitrary integration constant. The components of the spinor field can be 
obtained from (|2.34l) setting = and Q = (m — n) /2c. As far as F = F(J) is concerned, the 
volume scale x can be obtained from the corresponding equations setting spinor mass M = 0, while 
the for the components of the spinor field we need to set Q — (m — n)/2c and 4> = in (12.341) [In 
this case from (12.271) one finds P = Dq/x, with Do being some arbitrary integration constant]. 
Beside this, we assume c = c(x). The Eq. (12.491) then can be written as 

X = J?(x,q), (3.2) 

where we define 

J?(t, q) = 2(m 2 - mn + n 2 ) 4r + — \M + 3> + + 3At. (3.3) 

C L 2 L T' 

Here q is a set of problem parameters, namely, q = {m, n, M, X, T], A}. Equation (13.21 ) admits 
the following first integral, 

i=y/2[E-U(T,q)], (3.4) 

with the potential 

xdx 



U(x,q) = -{A(m 2 -mn + n 2 ) J '^ + 3k[Mx + J @dx+ T 1 ^] +3At 2 }. (3. 



5) 



From a mechanical point of view, Eq. (13.21) can be interpreted as an equation of motion of a single 
particle with unit mass under the force J?(x,q). In (|3.4I) E is the integration constant which can be 
treated as an energy level, and (x,q) is the potential of the force JP(q\,x). We solve Eq. (13.21) 
numerically using Runge-Kutta method. The initial value of x is taken to be a reasonably small 
one, while the corresponding first derivative x is evaluated from (13.41) for a given E. 
In what follows we solve Eq. (13.21 ) for (i) c = x and (ii) c = y/x. 



A. Case with c = x 

Let us consider this case in details for different types of spinor analyze field nonlinearity. 

1. Spinor field with power law induced nonlinearity 

Here we consider the case when the spinor field nonlinearity is given by a power law of S. In 
doing so we set F = S^, with rj being the power of nonlinearity. The right hand side of Eq. (13.21) 
in this case has the form 



, . m — mn + n 3 k 

&{x,q) = 2 + — 

X 2 



MC + , „ + 



+ 3At, (3.6) 



(Ti+ACj) 2 t? 
with C s = C 2 c Cq /2. For the potential in this case we have 

% (T, q) = -{4(m 2 - mn + n 2 )\n x + 3k[MC X - AC s /(t j? + XC% ) + CqT 1 ^] + 3At 2 }. (3.7) 
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We solve Eq. (13.21) with the right hand side given by (|3.6I) . For simplicity further we set ?c = 1, 
Co = 1, C s = 1 and £q = 1. For numerical calculations we used the following values for problem 
parameters: m = 2, n = 1 ,M = \,X = 0.1, £ = 1/3,7] = 3. Note that r\ = — 3 gives almost the same 
result as in case of r\ = 3. For cosmological constants we used A = 0, + 1, — 1, respectively. The 
energy level E is taken to be zero (E = 0). It should be noted that in the present case the potential 
possesses an infinitely high barrier at x = 0, it means in the case at hand x is always positive, that 
is we have singularity-free solution. Note that for a given value of E the minimum value of x 
(Tmin) should be greater or equal to the value of x at the point of intersection of E and % . Here 
we set E = 0, so T m i n > Tj nt : (T mt ) = 0. For simplicity here the initial value of x is taken to 
be unity. Moreover, in case of A < which is responsible for additional gravitational energy the 
value of x is bound from above as well, i.e., in this case the value of x lies between the two points 
of intersection of E and % . 




1 2 3 4 ^5 



FIG. 1: View of the potential given by (13.71) FIG. 2: Evolution of t for a negative A. 

for different values of A. 

In Fig. Q] we plot the potential corresponding to (13.71) . In Fig. [2] the evolution of x in case 
of a negative A is illustrated. As one sees, in this case the model allows an oscillatory mode 
of expansion. In Figs. [3] and 0] behavior of x for a nonnegative A is presented. As one sees, 
introduction of a positive A, which is often used to model the dark energy, results in the rapid 
growth of the universe. It is evident from the Eqs. (13.21) and (13.61) . the inhomogeneity plays 
crucial role at small x. In case of expanding universe inhomogeneity becomes notable only at the 
early stage of evolution, while for a oscillatory mode it is notable at x = T mm . Apparently from the 
expression (13.71) varying the value of inhomogeneity parameters m and n the height of the potential 
barrier in the vicinity of x = can be manipulated, though it will be infinitely high at x = 0. 

2. Spinor field with trigonometric induced nonlinecirity 

Here we consider the case when the spinor field nonlinearity is given by some trigonometric 
functions. Here we chose the following two cases with (i) F = s'm(S) and (i) F = exp(S). 
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FIG. 3: Evolution of x for a positive A. 



FIG. 4: Evolution of x for a trivial A term. 



If F = sin (5), we have 

_ , . „ m 2 — mn + n 2 3 

&(<C,q =2 + ^ 

T 2 



Acos 1/t l-£ 

^+^7 , , , H jr— ■ 

T 2 1+Asnl/T 2 



+ 3At, (3.8) 



and 



<%{x,q) = -{4(m 2 -mn + n 2 )lnT + 3fc[MT-A/2(l+Asin(l/T)) + T 1 "^]+3AT 2 }. (3.9) 
For F = exp(S) one finds 



. . „ m 2 — mn + n 2 3 
^ T, 9 =2 + x 

T 2 



Aexp(l/r) 1-C 



+ • 



T 2 (l + Aexp(l/faw)) 2 x^ 



+ 3At, (3.10) 



and 



<2f(T,$) = -{4(m 2 -mn + n 2 )lnT + 3?c[MT-A/2(l+Aexp(l/T)) + T 1 " ? ]+3AT 2 }. (3.11) 

Both these cases were solved numerically. The overall behavior of the potential % and x is 
almost the same as in case of a power law nonlinear term illustrated in the Figs. CE|-|4]). 



and 



B. Case with c = \fx 

In this case we find the following picture of the force J^"(t, q) and potential &(x, q): 

3 Atj ^ 1 — C 

^(x, q )=2(m 2 -mn + n 2 ) + -[M+ ^' + ^ 2 +— f-], 

ty{x,q) = -{4(m 2 -mn + n 2 )x + 3[Mx-X/(X + X 11 ) + X 1 -^}}. 



(3.12) 
(3.13) 
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Unlike the previous case now x can be trivial as well, thus giving rise to a spacetime singularity. 
Given the value of energy level E in case of a negative A the solutions may be both oscillatory and 
nonperiodic. 




1 2 3 4 5 6^7 



FIG. 5: View of the potential % given by (I3.13I ) FIG. 6: Evolution of x for a negative A with dif- 

fer different values of A. ferent values of E. 

Here as in previous case we set the following values for problem parameters: m = 2, n = 1 
M = \,X = 0.1, £ = 1/3,7] = 3. For cosmological constants we used A = 0, +1, — 1, respectively. 
For energy level E we set E = and E = — 1 . The evolution of x in case of A > is same as 
in previous case, while for A < the evolution of x is nonperiodic for E > and oscillatory for 
E < 0. In Fig. [5]we plot the potential corresponding to (13.131) . Contrary to the previous case it 
does not possess the infinitely high potential barrier at x = 0, which means in this case x can be 
trivial as well. In Fig. [6] we illustrated the evolution of x for a negative A with two different values 
of E. Given the concrete value of E the solution is either oscillatory or non-periodic. As far as 
inhomogeneity is concerned, in this case the part related to inhomogeneity can be added to the 
mass term and it plays the same role as the spinor mass in the evolution of the universe. 

IV. CONCLUSION 

A self-consistent system of interacting spinor and scalar fields within the framework of Bianchi 
type- VI (BVI) is studied in presence of a cosmological constant. Exact solutions of the spinor, 
scalar and gravitational field equations are obtained for some special choice of the spinor field 
nonlinearity. It is shown that introduction of a positive A which is often used to model the dark 
energy results in a rapid growth of the universe, while a negative A gives rise to an oscillatory or 
non-periodic mode of expansion. If the metric functions a and b are taken to be inverse to each 
other (ab = 1), we have a singularity free universe independent of the sign of the A term. 
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